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F-Term Hybrid Inflation Followed by a Peccei-Quinn Phase Transition 

G. Lazarides 1 ^ and C. Pallis 2 'Q 

Physics Division, School of Technology, Aristotle University of Thessaloniki, Thessaloniki 54124, GREECE 
2 Department of Physics, University of Cyprus, P.O. Box 20537, Nicosia 1678, CYPRUS 

We consider a cosmological set-up, based on renormalizable superpotential terms, in which a 
superheavy scale F-term hybrid inflation is followed by a Peccei-Quinn phase transition, resolving 
the strong CP and [s, problems of the minimal supersymmetric standard model. We show that 
the field which triggers the Peccei-Quinn phase transition can remain after inflation well above the 
Peccei-Quinn scale thanks to (i) its participation in the supergravity and logarithmic corrections 
during the inflationary stage and (ii) the high reheat temperature after the same period. As a 
consequence, its presence influences drastically the inflationary dynamics and the universe suffers 
a second period of reheating after the Peccei-Quinn phase transition. Confronting our inflationary 
predictions with the current observational data, we find that, for about the central value of the 
spectral index, the grand unification scale can be identified with its supersymmetric value for the 
relevant coupling constant re ~ 0.002 and, more or less, natural values, ±(0.01 — 0.1), for the 
remaining parameters. On the other hand, the final reheat temeperature after the Peccei-Quinn 
phase transition turns out to be low enough so as the gravitino problem is avoided. 

PACS numbers: 98.80.Cq 



I. INTRODUCTION 

One of the most natural and well-motivated inflation- 
ary model is the supersymmetric (SUSY) F-term hybrid 
inflation (FHI) [l], Q. It is realized at (or close to) 
the SUSY Grand Unified Theory (GUT) scale M GVT ~ 
2.86 x 10 16 GeV and can be easily linked to exten- 
sions [3[ of the Minimal Supersymmetric Standard Model 
(MSSM) which provide solutions to a number of prob- 
lems of MSSM. Namely, the ^-problem of MSSM can be 
solved via a direct coupling of the inflaton to Higgs super- 
fields [4j] or via a Peccei-Quinn (PQ) symmetry [5] which 
also solves the strong CP problem 6] . Also baryon num- 
ber conservation can be an automatic consequence [4J of 
a R symmetry and the baryon asymmetry of the universe 
can be generated via leptogenesis which takes place [7| 
through the out-of-equilibrium decay of the inflaton's de- 
cay products. 

The aforementioned resolution of the /i-problcm of 
MSSM via a PQ symmetry can be achieved Q by con- 
sidering non-renormalizable superpotential terms involv- 
ing additional singlets which develop vacuum expectation 
values (VEVs) of the order of the PQ symmetry break- 
ing scale. The emerging potential of these singlets has 
a local minimum which is separated from the global PQ 
minimum by a sizable potential barrier preventing a suc- 
cessful transition from the trivial to the PQ vacuum. As 
one can show [8|], by considering the one-loop tempera- 
ture corrections Q to the scalar potential, this situation 
persists at all cosmic temperatures after reheating. It is, 
thus, obligatory to assume that, after the termination of 
FHI, the system emerges with the appropriate combina- 
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tion of initial conditions so that it is led [10J to the PQ 
vacuum. Stated differently, the PQ symmetry is to be 
broken during or before FHI. As a consequence, tight up- 
per bounds on the inflationary scale have to be imposed 
so that the isocurvature fluctuations 11( of the axion are 
consistent with the observational bounds - this restric- 
tion, though, can be alleviated when cosmic strings are 
formed pjf . 

This latter complication can be avoided, if the PQ 
phase transition (PQPT) takes place after the end of 
FHI. This possibility can be realized adopting [13[ only 
renormalizable superpotential terms similar to those 
which lead to FHI. The fi parameter of MSSM can also 
be generated from the PQ scale as in Ref. [5]. On the 
other hand, this scheme may lead [141 ] to disastrous do- 
main walls which, however, can be avoided [l5l.ll6j by in- 
troducing extra matter superfields without jeopardizing 
the unification of the MSSM gauge coupling constants. 
Moreover, the structure of the superpotential obliges us 
to consider extra couplings between the inflaton field of 
FHI and the inflaton-like field of PQPT in the Kalher 
potential. As a bonus, these couplings assist us to recon- 
cile the results on the (scalar) spectral index n s of FHI 
and the recent seven-year results [ItJ from the Wilkinson 
microwave anisotropy probe (WMAP7) satellite. 

Indeed, it is well-known that the realization of FHI 
within minimal Supergravity (SUGRA) leads to n s which 
is just marginally consistent with the fitting of the 
WMAP7 data by the standard power-law cosmological 
model with cold dark matter and a cosmological constant 
(ACDM). O ne p ossible resolution (for other proposals, 
see Refs. HHllf) of this problem is [H, Hi] the addition 
to the Kahler potential of a non-minimal quatric term of 
the inflaton field with a convenient choice of its sign. As 
a consequence, a negative mass term for the inflaton is 
generated. In the largest part of the parameter space, 
the inflationary potential acquires a local maximum and 



minimum. Then, FHI of the hilltop type [25] can occur 
as the inflaton rolls from this maximum down to smaller 
values. Therefore, n s can become consistent with data, 
but only at the cost of an extra indispensable mild tuning 
|23| of the initial conditions. Another possible complica- 
tion is that the system gets trapped near the minimum 
of the inflationary potential and, consequently, no FHI 
takes place. 

Within our proposal, a negative mass term for the in- 
flaton is generated too, which however depends on the 
various coefficients (of the order 0.01 — 0.1) involved in 
the Kahler potential and not exclusively on the quatric 
term of the inflaton field. As a consequence, sufficiently 
low n s can be achieved with cither sign of the terms in the 
Kahler potential and the inflationary potential remains 
monotonic during the period of FHI. Therefore, compli- 
cations due to the appearance of maxima and minima 
along the inflationary path can be eluded. Another by- 
product of our proposal is that the reheat temperature 
after PQPT, calculated consistently with the model in- 
teractions, turns out to be low enough so that the grav- 
itino (G) constraint [2a. |27| and the potential problem of 
topological defects |28| of FHI can be significantly relaxed 
since an entropy release takes place after PQPT. On the 
other hand, the reheat temperature is high enough so 
that non-perturbative electroweak sphalerons are opera- 
tive and, consequently, thermal [29j or non-thermal [30j 
leptogenesis can, in principle, work. 

Below, we present the basic ingredients of our model 
(Sec. (Til and describe the inflationary potential and dy- 
namics (Sees. Mil and HVj) . We then exhibit the con- 
straints imposed on our cosmological set-up (Sec.|V|. We 
end up with our numerical results (Sec. IVI|) and our con- 
clusions (Sec. IVIIj) . 



II. MODEL DESCRIPTION 
A. The general set-up 

In order to explore our scenario, we adopt the 
left-right symmetric gauge group Glr = SU(3) c x 
SU(2) L x SU(2) R x U(l) B _ L . In our scheme, Glr can 
be broken down to the standard model gauge group Gsm 
at a scale close to the SUSY GUT scale Mgut through 
the VEVs acquired by a conjugate pair of SU(2)r dou- 
blet left-handed Higgs superfields, $ and $, with B — L = 
— 1, 1 respectively. As a consequence, no cosmic strings 
are produced [28|, [31] in this realization of standard FHI 
and, therefore, we are not obliged to impose extra restric- 
tions on the parameters - as e.g. in Ref. |32j |. 

The model possesses also three global U(l) sym- 
metries. Namely, a (color) anomalous PQ symmetry 
U(1) P q, an anomalous R symmetry U(1)_r, and the 
baryon number symmetry U(l) s . Note that global con- 
tinuous symmetries can effectively arise [33j j from the rich 
discrete symmetry groups encountered in many compact- 
ified string theories (see e.g. Ref. [3J]). The PQ sym- 



TABLE I: Superfield Content of the Model 



Superfields Representations 
under Glr 



Global Symmetries 
R PQ B D 



Matter Fields 



li 


(1,2,1,-1) 





-2 










(1,1,2,1) 


2 











Qi 


(3,2,1,1/3) 


1 


-1 


1/3 





it 


(3,1,2,-1/3) 


1 


-1 


-1/3 






Higgs Fields 

S (1,1,1,0) 4 

$ (1,1,2,-1) 

J> (1,1,2,1) 



P 

Q 

Q 



(1,1,1,0) 
(1,1,1,0) 
(1,1,1,0) 



4 
-2 
2 



(1,2,2,0) 



Extra Matter Fields 



A, (3,1,1,2/3) 2 1 0-1 

D a (3,1,1,-2/3) 2 10 1 

H a (1,2,2,0) 2 10 



metry U(1) P q can be spontaneously broken at the PQ 
breaking scale f a ~ (10 10 — 10 12 ) GeV (which coincides 
with the axion decay constant - for a review see Ref. [35]) 
via the VEVs acquired by two Glr singlet left-handed 
superfields Q and Q. 

The part of the superpotential which is relevant for the 
breaking of the group Glr X U(1) P q is 

W = kS (** - M 2 ) + K a P(QQ - fl/A) + XSQQ, (l) 

where S, P are Glr singlet left-handed superfields which 
trigger the breaking of Glr and U(1) P q respectively and 
the parameters k, K a , M ~ Mgut, and f a are made pos- 
itive by field redefinitions. For simplicity, we restrict our 
analysis to real A's. The superfield P can be regarded 
as the linear combination of the Glr singlets with PQ 
and R charge equal to and 4 respectively that does not 
couple to $$ - cf. Ref. [36|. In this basis, the most gen- 
eral Kahler potential of our model includes interference 
terms of S and P even at the quadratic level - contrary 
to the choice opted in Ref. |37}. Namely, we adopt the 
following Kahler potential 

K = \S\ 2 + \P\ 2 + a(SP* + S*P) 



b \sl 

4mp 

e|5| 2 - 



JPj* 
'4mp 



,\S\ 2 \P\ 



2m 



f\P\ 2 

tl—L (SP* + S*P) 



4m 2 



(sp*y + (s*py 
m 2 + m 2 + \Q\ 2 + \Q\ 2 

k($$ + $*$*) + l{QQ + Q*Q* 



, (2) 



where all the coefficients a, b, c, d, e, /, g, k, and I are 
taken, for simplicity, real although some of them can be, 
in general, complex. The ellipsis represents higher or- 
der terms involving the waterfall fields ($, $, Q, and Q) 
and S and P. We can neglect these terms since they have 
negligible impact on the SUSY vacuum and are irrelevant 
along the inflationary path - see below. 

The usual superpotential terms of MSSM can be de- 
rived from the following superpotential: 



2 1,2 



W n 



A, 



Q 2 h 
2tot 



*j?*^ 

Vvij 1- yiijkhlj + y q i 3 qihq c 3 (3) 



with mp ~ 2.44 x 10 18 GeV being the reduced Planck 
scale. Here, the ith generation SU(2)l doublet left- 
handed quark and lepton superfields are denoted by qi 
and li respectively, whereas the SU(2)r doublet anti- 
quark and antilepton superfields by qf and l\ respec- 
tively. The electroweak Higgs superfields are contained 
in a SU(2) L x SU(2) R bidoublet Higgs superfield h. The 
first term in the right-hand side (RHS) of Eq. © gen- 
erates the n term of MSSM via the PQ breaking scale 
(see below), while the second term generates intermedi- 
ate scale masses for the right-handed neutrinos and, thus, 
seesaw masses |2| for the light neutrinos. 

The representations under Glr and the charges under 
the global symmetries of the various matter and Higgs 
superfields contained in this model are presented in Ta- 
ble HI which also contains the extra matter superfields 
required for evading the domain wall problem associated 
with the PQPT together with a new imposed global sym- 
metry U(1) D - see Sec. uTCl 



B. The cosmological scenario 

The F-term SUGRA scalar potential obtained from W 
in Eq. (TTJ) and K in Eq. ^ can be found by applying the 
well-known formula (see e.g. Ref. [l|) 



Vsugra = e 



k i^{f:,{k^)- 



F* 



,|Wf 



(4) 



where Fi 



W A 



K t iW/rrip and a subscript ,i [,'< 



denotes derivation with respect to (w.r.t.) the complex 
scalar field i [i*]. In the limit where mp tends to in- 
finity, we can obtain the SUSY limit, Vp, of Vsugra in 
Eq. ((4]) which turns out to be 

U F ~ | K ($$-Af 2 ) +XQQ\ 2 /{1-a 2 ) 
-J 2 /4| 2 /(l-a 2 ) 



m 2 



i$n 



K 2 \S\ 2 y 

\XS + Ka P\ 2 (\Q\ 2 + \Q\ 2 ) 

^a (Q*Q* - fa/4) [K (** - M 2 ) 

XQQ]/(1 — a 2 ) + complex conjugate , (5) 



where the complex scalar fields which belong to the SM 
singlet components of the superfields are denoted by the 
same symbol and we take into account that M ^> f a . 
Note that although the mixing terms in the first line of 
Eq. ([2]) have an impact on Eq. ([5]), the corresponding 
terms in the last line of this equation do not contribute 
since they are holomorphic or anti-holomorphic in the 
field variables. Also the terms in the ellipsis in Eq. @ 
have no appreciable contribution to the RHS of Eq. ([5]) 
since they are suppressed by powers of mp . 

The D-term contribution vanishes along the direction 



1*1 = 1*1. 



(6) 



From the potential in Eq. ([5]) and taking into account 
that M ^ f a , we find that the SUSY vacuum lies at 

(S) = 0, |<*)| = |<*)| ~ A/, (P) = 0, and \{<t> Q )\ = f a , 

where we have introduced the canonically normalized 
scalar field 4>q = 2Q = 2Q. Note that, since the sum 
of the arguments of (Q), (Q) must be 0, Q and Q can 
be brought to the real axis by an appropriate PQ trans- 
formation. As a consequence of Eqs. Q, W leads to a 
spontaneous breaking of Glr and Upq(1). The same su- 
perpotential W gives also rise to a stage of FHI and a 
PQPT, which is highlighted in the following. 

The potential in Eq. © possesses a D- and F-flat 
direction at 



$ = $ = 0andQ = Q = 
with a constant potential energy density 

K 2 M 4 



Vhio 



1-a 2 



(since M > /„). 



(8) 



(9) 



The direction in Eq. ([8]) can be used as the inflationary 
trajectory since it corresponds to a classically flat valley 
of minima for 



(a) \S\ > -^L= and (b) \a a \ > J " {X ^ a) M. 



VT 



\-& z 



(10) 
where we have defined o~ a = AS 1 + n a P. 

For relatively large A's, the reheating temperature af- 
ter FHI turns out to be rather high, leading to troubles 
with the G abundance [26|, [27| - as usually in SUSY cos- 
mological models. However, we can verify that a D- and 
F-flat direction appears at 

S = 0, * = * = M, and Q = Q = (11) 

with potential energy density 

^pqo = ^/ a 4 /16. (12) 

Since Vpqo <C Uhio: ^pqo can temporally dominate over 
radiation after the end of FHI and, if \P\ > f a , drive a 
PQPT - see Sec. HVi A subsequent second episode of 



TABLE II: The mass spectrum of the model along the infla- 
tionary trajectory of Eqs. ([5| and (|10l) . 



Superfields 
of origin 


Fields 


Mass Squared 


Bosons 


Q,Q 


4 complex scalars 
2 complex scalars 


|2 , K(A-are a )M 2 


Fermions 


Q,Q 


2 Dirac spinors 
1 Dirac spinors 


k 2 |S| 2 

Kl 2 



reheating can dilute sufficiently the unwanted G concen- 
tration - see Appendix |XJ At the end of PQPT, the 
fields P, Q, and Q acquire their VEVs in Eq. Q. The 
\i term of the MSSM is generated via the first term in 



the RHS of Eq. © with 



Ml 



■V (Q) / to p> which is of 



the right magnitude if \{Q)\ = f a /2 ~5x 10 11 GeV and 
A M ~ (0.001-0.01). 

The cosmological scenario above can be attained if we 
ensure that, at the end of FHI, the fields $ and $ acquire 
their VEVs in Eq. (J7]), while the fields Q and Q remain 
equal to zero. To quantify this crucial requirement, we 
construct the mass spectrum of the theory along the in- 
flationary path of Eqs. ((5J and (|T0|) . We summarize our 
results in Table [TTJ We see that the mass matrices of 
the scalar components of the $, $ and Q, Q superfields 
develop a negative eigenvalue as |5| and \a a \ cross below 
their critical values - i.e. their lower bounds in Eq. (1101) . 
However, if the tachyonic instability of the $ — $ system 
occurs first, <! and <& start evolving towards their VEVs, 
whereas Q and Q continue to be confined to zero. 



C. Evading the domain- wall problem 

It should be mentioned that instanton and soft SUSY 
breaking effects explicitly break U(l)^ x U(1)pq to a dis- 
crete subgroup. Spontaneous breaking of this subgroup 
at the PQPT by (Q) and iff) can lead [14J to a disastrous 
domain wall production since this transition occurs after 
FHI. In order to avoid this problem, we must introduce 
some extra matter fields (see Table [TJ following Ref. [15| 
(see also Ref. [16| . where the set-up is similar to ours). 
Namely, we introduce n pairs of left-handed superfields 
D a and D a (a = 1, ...,n) which are SU(3) C antitriplets 
and triplets respectively with R charge equal to 2 and 
PQ charge equal to 1. However, these fields acquire in- 
termediate scale masses after the PQ breaking, which 
could prevent the unification of the MSSM gauge cou- 
pling constants. To restore gauge unification, we include 
an equal number of SU(2)l X SU(2)r bidoublet super- 
fields H a with PQ and R charges equal to those of D a 
and D a . In accordance with all the imposed symmetries, 
we can give intermediate scale masses to D a — D a and 



H a through the superpotential terms 

W dw = \ Da QD a D a + X Ha QH 2 a . 



(13) 



Here, we chose a basis in the D a — D a and H a space where 
the coupling constant matrices Au a and A# a are diagonal. 
Note that the full superpotential is invariant under a new 
global \](1) D symmetry - see Table U During FHI, we 
take D a — D a = H a = and, at the SUSY vacuum, we 
have (D a ) = (D a ) = (H a ) = 0. 

The number n can be determined as follows: The ex- 
plicitly unbroken subgroup of U(1).r x U(1)pq can be 
found, for every n, from the solutions of the system 

4a = (mod 2tt) and - 12a + 2(n- 6)/3 = (mod 2tt) , 

(14) 
where a and /3 are the phases of a U(l)ij and U(1)pq 
rotation respectively. Here we took into account that the 
R charge of W and, thus, of all the soft SUSY breaking 
term is 4 and that the sum of the R [PQ] charges of the 
SU(3) C triplets and antitriplets is —12 [2(n — 6)]. We can 
induce, therefore, that the explicitly unbroken subgroup 



is Za x Z. 



2(n-6)' 



It is then important to ensure that this 



subgroup is not spontaneously broken by (Q) and (Q), 



i.e. 



2/3 = (mod 2tt) 



(15) 



since otherwise cosmologically disastrous domain walls 
will be produced at PQPT. This requirement implies that 
n must be 5 or 7. In both these cases, the subgroup of 
U(l)ii x U(1)pq left unbroken by instantons and SUSY 
breaking coincides with the one left unbroken by (Q) and 
iff) and is & Z4 x Z2- 

It is easy to find the contribution of Z) a , D a , and H a to 
the coefficients b\ , 62 , and 63 controlling [38| the one loop 
evolution of the three gauge coupling constants (71,(72, 
and 53 within MSSM. It is then straightforward to prove 
that, if we assign B — L = 2/3 and —2/3 to D a and D a 
resp ectively, the quantities 62 — b\ and 63 — £>2 (which are 
[38| crucial for the unification of 31,32, and 33) remain 
unaltered. Therefore, the inclusion of the extra matter 
superfields does not disturb the gauge unification at one 
loop. 

Recapitulating this section, let us present the total su- 
perpotential of our model which is 



w tot = w + w n 



w dv 



(16) 



where W, W m and Wd w are given in Eq. ([T]), (J3J) , and 
(fP3")) respectively. 



III. THE INFLATIONARY EPOCH 

A. Structure of the inflationary potential 

The potential which can drive the inflationary stage of 
our set-up has the following general form 



Vhi — Vhio + Vkls + VhIc, 



(17) 



where Vnis and Vkic represent, respectively, SUGRA and 
one-loop radiative corrections to the inflationary poten- 
tial, calculated in Sec. MI A II and IIII A 21 Let us note, 
in passing, that the most important contribution [39] to 
Vhi from the soft SUSY breaking terms is expected to 
start playing an important role for rather small k's and 
so, it remains negligibly small in our set-up due to the 



J 



large k's used - see Sec. IVT1 



1. Supergravity corrections 

Expanding the F-term SUGRA potential in Eq. (gj) 
along the inflationary trajectory - see Eq. © - in powers 
of 1/mp for M ^ f a , we obtain the following expression: 



Uhis 



Vt 



HIO 



(1 - a 2 )mp 



^i|S| 2 + A 12 (S*P + PS*) + A 2 \P\ 



Vj 



HIO 



4(1 - a 2 ) 2 m 4 
B 3 \S\ 2 \P\ 2 + (B 4 \S\ 2 + B 5 \P\ 2 ) (S*P + PS*) + B 6 ((S*P) 2 + (P*S) 2 ) 



B^Sf + B 2 \P\ 4 



where 



(18) 



Ax = 2ae - a 4 - a 2 (d - 1) - 6, A 2 = 1 - d - a(a + ac - 2/), A 12 = a(l + d - a(a + /) + g) - e, (19) 

Si = 2 + a 4 (2 + Ad - 36) + 6(46 - 7) + 4a 5 e - 8a 3 (2 + d)e - 4a(2d + 46 - 3)e + 4e 2 

+2a 2 (56 + 2d(d+26- l) + 6e 2 -2), (20) 

B 2 = 2 + c + a 4 (5c + 2) + 4(d - l)d - 8a 3 (c + 1)/ - 8a(c + 2d - 1)/ + 4/ 2 + 2a 2 ((2c + l)(c + 2d - 2) + 6/ 2 ),(21) 
P 3 = 4[1 + a 6 + (d - 2)d + b{2d - 1) + 4a 5 / + 2e(e + /) - 2a(ce + 2(6 - 1)/ + (2e + f)(2d + g)) + g 2 

+a 4 (c - 2(d + 5) - l) - 2a 3 (ce + /(4 + 2rf + gi)) + a 2 (6(1 + 2c) + 4rf + c(2d - 1) + 2/(5e + /) 

+2g + (d + g)(3d + g)-l)], (22) 

S 4 = 2[2a 4 (/ - 3e) + 2e(26 + d + g - 2) + a 5 (2 + 4d + .g) + 2a 3 (6 - d(4 + d + g) - 2(1 + ef + g)) 

+2a 2 ((26 + 2d - 1)/ + e(5 + bd + 3g)) - a(4e(2e + /) - 3g + 2d{d + g - 2) + 2b(2d + 2g + 1) - 2)] , (23) 
B 5 = 2 [2 ((2d - l)e + /(d + g)) - 2a 4 / - 2a 3 (2/ 2 + g + c{d + g + 2) + 2) + 2a 2 (e + 2ce + /(2c +5d + 3g+ I)) 

+a 5 (2 + 3c) + ae(l - 2d - 2 5 ) + 2a(l - 2/(2e + /) + g - 2d(d + g))] , (24) 



and 



B 6 = 2a 6 + 4e 2 + 6a 5 / + (4d-l)g- a 4 (4 + 4d + 5g) + 4a 3 (e - /(2 + d + g)) - 2a(/(2d + 2.g - 1) 
+e(2 + 4d + 4.g)) + 2a 2 (2/(2e + /) + 3.g + 2(d + d 2 + gd + g 2 ) + l) . 

I 



(25) 



In this expansion, we neglect terms proportional to first term of Vnis in Eq. (fl~8| reads 



y/VmoVpQo/rrip, which are numerically suppressed. To 
specify the canonically normalized scalar field, which 
could play the role of inflaton, we have firstly to bring 
into a canonical form, through a non-unitary transfor- 
mation, the quadratic part K$p of the Kahler potential 
in Eq. ([2]) which involves S and P. Namely, we find 



K 



SP 



\Sf 



LPI 2 - 



-a{SP*+S*P) 



(26) 



where s± = A± (S ± P) fy/2 with A± = y/l ± a. Solving 
w.r.t. S and P, we get 



S 



W-4 

V2 Uh 



S- 

A~l 



and P 



1 

71 



A~l 



(27) 

Introducing the real and imaginary components of the 
fields s± from the relations s± = {s%± + is 2 ±)/\/2, the 



Uhio 
2m 2 



(t s - *+)*#[',») + (i ^4 



where the matrices M 2 and M 2 are found to be 

Ai+A2-2A 12 Ai-A 2 

Ml 



Ml 



(l_ a )2( 1+a ) 

A t -A 2 

(1-a 2 ) 3 / 2 



(l_ a 2)3/2 



(28) 



(29) 



(l+a)2(l-a) 

To identify the combination of s± which can play the role 
of inflaton, we have to diagonalize the matrices M 2 and 
M 2 ■ This can be realized via an orthogonal matrix U as 
follows (t stands for the transpose of a matrix): 



UM^U T = UM 2 U T = diag (m 2 + ,m 2 _ 



where 



m\. 



D 1 ±^D 2 



(30) 



(31) 



with 

and 

Also 

with 

U± = 



Di = Ai- 2aA 12 + A 2 



D 2 = D 2 + 4(1 - a 2 ){A 2 2 - AtA 2 ) 



U 



f 1/N+ 1/N_ ) 
{U+/N+ V-/N-) 



T +/N+ U- 
aAx - 2A 12 + aA 2 ± VAi' 



(32) 

(33) 
(34) 



VT^^{A 1 - A 2 ) 



and N± 



i + ul. 



(35) 



Embedding unity (1 = UU T = U T U) on the left and the 
right of M 2 and M 2 in Eq. ([28]) . this expression can be 
brought into the form 



^HIO 



[m\\(j+\ 



2 i 2 I |2 

+ m_ ct_ 



(36) 



where the real and imaginary components of the complex 
fields cr_|_ and tr_, defined as <r± = (er^ + ia 2 ±)/\/2, are 
given by 

o-i- = Si-/N + + Si+/N_ and (lo2), (37) 

o-i + = Sl -U+/N + + si+i7_/iV_ and (1 <-> 2). (38) 

Performing an appropriate R transformation, we can ro- 
tate <7_ to the real axis (setting a 2 - =0). To simplify 
the notation, we rename the remaining fields as follows: 
cti_ = er, (Ti + = s and (72+ = <?• We then can solve 
Eqs. ([37| and (|38]l w.r.t. the components of s±, i.e., 



Sl- 



-<V+ — t:— and S2- 



JV+g 



«i+ = N - 



U+-U- 
s — U+o~ 

U--U+ 



and s 2 + = 



U+-U- 

N_q 
U- -C/ 4 



(39) 
(40) 



Substituting the expressions above into Eq. (f2~Tf and then 
into Eq. (|18p . we can derive Vhi s as a function of a, s, and 
q. As can be explicitly verified, K$p remains canonical 
w.r.t the latter fields - i.e. K S p = (a 2 + s 2 + q 2 ) /2 - 
since the U transformation, which connects <7, s, and q 
with si± and s 2 ±, is orthogonal. 

From the expressions above, we easily conclude that, 
if we set b = c = d = e = f = g = with a ^ 0, we take 



-4i 



A 



12 



-4, 



1 



J 2 'l 



f 



(a 2 ,a,l), 



(41) 



which results to m_ — and to+ = 2 via Eq. pip . In 
other words, in this case, a remains identically mass- 
less, whereas s and q acquire effective masses equal to 
\/3/(l - a 2 )ff H io with iJ H io = VVmo/V^mp. There- 
fore, we are obliged to invoke non-zero coefficients of the 
higher order terms of the Kahler potential considered in 



Eq. ([2]), in order to generate a negative to 2 ,, and to re- 
duce, thereby, n s to an observationally acceptable level. 
Needless to say that, for a — b — c = d — e = f = g — 0, 
the well-known results (see e.g. Refs [4fJ, |4J|) in the 
context of minimal SUGRA can be obviously recovered 
- note that, in this case, B\ = B 2 — B%/2 = 2 and 
B4 = B$ = Bq = 0. 



2. Radiative corrections 

The constant tree-level potential energy density Vhio 
causes SUSY breaking leading Q to the generation of 
one-loop radiative corrections, which provide a logarith- 
mic slope along the inflationary path necessary for driv- 
ing the system towards the vacua. Inserting the spectrum 
shown in Table |TT] in the well-known Coleman- Weinberg 
formula |43|, we find (compare with Ref. [36|]) that the 
one-loop radiative correction to Vki is 



Vt 



Hie 



V K + V X , 



(42) 



where 



with x = (1 - a 2 )\S\ 2 /M 2 and 



(A-aK a ) 2 U H io /„, k(X - aK a )x a M 2 

V * = o^ 9/, -^- I 2 m ~ ^-T5 1" Ac {X a ) 



32tt 2 (1 - a 2 ) 



(1 - a 2 )A 



2U2 



with x a = (1 — a 2 )|cr a | 2 /K(A — aK a )M 2 . In the relations 
above, we have taken into account that the dimensional- 
ity of the representations to which <t and $ [Q and Q] 
belong is 2 [1] - see Table I - and we have defined 

f rc (y) = (y+l) 2 ]n(l + l/y) + (y-l) 2 ln(l-l/y). (45) 

Let us note that, for x ^> 1 and x a 3> 1, V K [V\] can be 
well approximated by replacing f rc (y) — 3 in the RHS of 
Eq. Q [Eq. dH])] - see Sec. lnTBl Although rather large 
k's, k q 's, and A's are used in our work, renormalization 
group effects [43j remain negligible and so our results are 
independent from the renormalization scale A. 

Substituting Eqs. ^ and flm into Eq. $27$ and then 
into Eqs. (|4*3")l and (pH]) , we can derive Vnic as a function 
of cr, s, and q. 



B. The inflationary dynamics 

The evolution of the various fields involved in our 
scheme during FHI is governed by their equations of mo- 
tion (e.o.m.) 



/ + 3Hf + Vhi,/ = with / = a, s, and q. 



(46) 



where the dot denotes derivation w.r.t. the cosmic time 
t and H is the Hubble parameter. The solution of the 



system in Eq. (|46j) can be facilitated if we use as inde- 
pendent variable the number of e-foldings N denned by 

TV = In (R/Rwi) => N = H and H = H'H (47) 

with R{t) being the scale factor of the universe and Ran 
its value at the commencement of FHI. Here the prime 
denotes derivation w.r.t. N. Converting the time deriva- 
tives to derivatives w.r.t. N, Eq. (|46p becomes 

H 2 f" + {3H + H')Hf + Vmj = (48) 

with / = c, s, or q. This system can be solved numeri- 
cally for the period of FHI by taking 



Since the general solution of the corresponding homoge- 
nous differential equation rapidly decreases to zero, as the 
field q, the solution of Eq. (|48| with / = s, is dominated 
by the following particular solution 



2k 2 Ci + (A - aK a ) 2 C 2 
8(1 - a 2 )n 2 m 2 + a ' 



/p. 



(54) 



which minimizes Vhi in the s-direction, as can be seen 
from Eq. (151| . In sharp contrast to the situation of 
Refs. [441446] . in our case s turns out to be just mildly, 
and not drastically, reduced w.r.t. a. 

In the slow-roll approximation, which is determined by 
the condition 



H = m 



2Vkio 



HI 



6mp - a' 2 - s' 2 - q' 2 '' 



(49) 



V = Vhi given by Eq. (fT71) . and imposing the following 
initial conditions (at N = 0): 

/(0) = fun and /'(0) = with f = a, s, or q, (50) 

where /hk is taken to be in the range (1.5 — 4.5) x 
10 17 GeV. We checked that our results are pretty stable 
against variation of /mi - see also Sec. I VI I 

Nevertheless, we can obtain a comprehensive and 
rather accurate approximation to the inflationary dy- 
namics if we put H ~ Hhiq = VVmo / VSrnp (H' = 0) 
and keep in Eq. (|48j) the most important terms in the 
expansion of Vhi,/- In particular, we can easily verify 
that Vhi,,} turns out to be proportional to q. As a conse- 
quence, Eq. (|48|) for / = q is a second order linear homo- 
geneous differential equation which admits an oscillatory 
solution with decreasing amplitude. Therefore, q rapidly 
decreases to zero and so it does not influence the dynam- 
ics of the other fields. On the contrary, the e.o.m. of 
s is non-homogenous as we can realize by inserting into 
Eq. P5|) for / = s the expression 



Vhi, s — Vhio 



2K 2 d + (A - an a ) 2 C 2 

8(1 - a 2 )TT 2 <T 



> (51) 



where the last two terms in the RHS of this equation have 
been derived by expanding the exact result for sj a — > 0. 
Here 

r _* A 2 JV S + n 2 a N P + 2\ Ka N SP 

D s X 2 D S + KiDp + 2\K a D S p 

where we have used the abbreviations: 

N\U- N 2 U+ N_N + (U- + U+) 



N s[P] - 



A 2 _ A\ 

r 2 7V Af2 1 



N'iU+ N%U- N'iU'i NlUl 

NsP — 75 75 — , D SP — 



A+A. 
NtUl 



A\ 



A 2 



A\ 



A 2 



D 



S[P] 



N+U- 
A- 



N-U+ 

~aT 



(53) 



max{e(cr),|?7(cr)|} < 1, 



(55) 



where 



(Vm,o 



2 V Vhi 



the e.o.m. of a in Eq. 
Vhi,<t or 



and r\ 



Vhi j0 



takes the form —3H 2 a' 



(56) 



2 mp(2« 2 + (A — &Ka) ) 

m_a + 



C 3 a 3 



8TT 2 (l-a 2 )a 



mp(l — a 2 ) 2 



where - keeping only the most important terms in the 
expansion of the second term of the RHS of Eq. (Q2 
we have 



(58) 



c , 


/d n2 i o_ 7~>2 


Ca - 16(U_-U + )*^ X " S ' "*" p 


+ B 3 D 2 SP + 2B 4 D SP Ds). 


Solving Eq. ([BT)l w.r.t. a, we get 


m P (l — a 2 ) 

o" = , — 

\/2^ 


- ml - D tan ( DN 


1 

— arctan — 

D 


V - m 2 p (l - a 2 ) 2 J J 



1/2 



(59) 



with (Thu being the initial value of a at the beginning of 
FHI and 



D= ^C 3 (2k 2 + (A - a Ka ) 2 )/27r 2 (l - a 2 ) 3 



(60) 



The end of FHI takes place at a = at , when the condi- 
tions of Eq. (P2J or Eq. d55j) are violated. For k < 0.005, 
one can consider the former case only since the slow-roll 
conditions are violated infinitesimally close to the rele- 
vant critical point of the inflationary trajectory. 

Soon afterwards, the inflaton system consisting of the 
two complex scalar fields S and (J$ -I- J$)/\/2 (where 
c>i> = <J — M and 5$ = $ — M) with mass mimf = V2kM 
settles into a phase of damped oscillations about the 
SUSY vacuum and decays to MSSM degrees of freedom 



8 



reheating the universe. The predominant decay channels 
of S and (<5i > +<5$)/v / 2 are to fermionic and bosonic Q, Q 
respectively via tree-level couplings derived from the last 
term in the RHS of Eq. ([!]) with a common decay width 



T^ \2 2 



(61) 



The corresponding reheat temperature is [471 ] given by 



Tlrh 



72 



5TT 2 girh* 



1/4 



y/Ti 



m P , 



(62) 



where gi Y h* = ff*(Tirh) counts the effective number of 
relativistic degrees of freedom at temperature T\Th- We 
find 5i rh » ~ 438.75 [g llh * ~ 513.75] for the MSSM spec- 
trum plus right handed neutrinos and the particle con- 
tent of the superfields P, Q, Q, D a , D a , and H a for n = 5 
[n = 7]. 

For relatively large A's (~ 0.05 — 0.1), we get Ti r h > 
VpQ . As a consequence, after the end of FHI, we ob- 
tain matter domination (MD) for T > Ti r h and radia- 
tion domination (RD) for V p q < T < Tirh. During MD 
the evolution of s can be found by solving its e.o.m. for 
N>Nw, i-e., inserting dHHEl 



H = H Hw e-^ N - N ^/ 2 , V = 



3H 2 



(63) 



into Eq. (J48J). Taking also into account that, during the 

-I /o 

MD epoch, R oc p OS c (where p sc is the energy density 
of the oscillating system), we can derive the value of s at 
the beginning of PQPT: 



SPQi 



( Plrh 

\Vkio 



1/4 



Shh with p llh = — girh*Ti Th (64) 



being the radiation energy density at temperature Ti r h 
and shh the value of s at the end of FHI, which can 
be approximated by inserting a = at into Eq. (|54[) . Note 
that during the subsequent RD era, s remains frozen (see 
Ref. [41| and footnote 1 in Ref. [4y|) and so the further 
reduction of s relative to a can be evaded - cf. Ref. 144| . 



IV. THE STAGE OF PQPT 

For T < Vfjio i the cosmological dynamics is governed 
by the second term of W in the RHS of Eq. (fJJ). The 
SUGRA corrections can be safely ignored since \P\ ~ 
/„ -C ttt-p- The relevant F-term scalar potential is 



potential - cf. Sec. IIII A 21 Indeed, the particle spec- 
trum there includes a Dirac spinor with mass K a |P| and 
2 complex scalars with mass squared k 2 (\P\ 2 ± /„/4) 
and, therefore, Vpq c can be written [42[ as follows: 



VpQc 



32tt 2 



21n- 



A 2 



+ /rc(»») 



where 



Ap = 



A_ AT_ + A+ N + 
A + A_{U.-U+) 



(66) 



(67) 



with x s = Aps/f a and P = A P s/2 - see Eqs. fl27|), (|39|) . 
and PDI . When \P\ < f a /2, one mass squared becomes 
negative and suggests a phase transition along the \P\- 
axis. Assuming gravity mediated soft SUSY breaking, 
the potential there has the form 



^PQ 



^PQC) + o m 



V 



^2V/pqo |a s |s + y PQc (68) 



for s > f a /Ap. Here, m s is the soft SUSY breaking mass 
of s and a s is the soft SUSY breaking tadpole 39]. In 
Eq. (|SB]h we take the phase of a s to be arg(a s ) = n which 
minimizes the Vpq for given s > 0. For reasonable val- 
ues of the parameters involved, Vpq does not give rise to 
inflation, mainly due to the contribution of Vpq c which 
spoils the 77-criterion - and has not been taken into ac- 
count in the analysis of Ref. |1J. This negative result 
remains even if we consider dissipative effects [50j due to 
the presence of the extra superpotential terms in Eq. (|13l) . 
Therefore, we are obliged to assume that after a negli- 
gible number of e-folds, the system consisting of the two 
complex scalar fields P and (SQ + SQ)/\/2 (where SQ = 
Q-f a /2 and SQ = Q-f a /2) with mass to pq = K a f a /\/2 
enters into an oscillatory phase about the PQ minimum 
and eventually decays, via the first non-renormalizable 
coupling in the RHS of Eq. (|3|) , to Higgses and Higgsinos 
respectively with a common decay width [3| 



r 2 = -L\ 2 (Ja- 

2-k M V 2m P 



m-PQ- 



(69) 



The corresponding reheat temperature is [47| calculated 
by 



T2rh — 



72 



57r 2 i g 2 rh* 



1/4 



V^ 



imp, 



(70) 



where ^2rh* = 228.75 is the energy density effective num- 
ber of degrees of freedom for the MSSM spectrum. 



V. OBSERVATIONAL CONSTRAINTS 



^pqf = n 2 a \QQ - f a /4\ 2 + n 2 a \P\ 2 (|Q| 2 + |Q| 2 ) . (65) 

Along the F-flat direction in Eq. (fTTj) . Upqf takes the 
constant value in Eq. (|12l) . which breaks SUSY creating 
a mass splitting in the supermultiples Q and Q and giving 
rise to one- loop radiative correction, Vpq c , to the relevant 



Under the assumptions that (i) the curvature perturba- 
tion generated by a is solely responsible for the observed 
curvature perturbation and (ii) the violation of Eq. (|10p 
occurs along the lines of Sec. Ill BL the parameters of our 
model can be restricted imposing the following require- 
ments: 



a. According to the inflationary paradigm, the 
horizon and flatness problems of the standard Big 
Bang (SBB) can be successfully resolved provided that 
the number of e- foldings, Ahi*, that the scale k* — 
0.002/Mpc suffers during FHI takes a certain value which 
depends on the details of the cosmological scenario. Em- 
ploying standard methods [20|, we can easily derive the 
required Ahi* at fc*: 



2 V 1/4 



Vr 



1/4 



5 ln iG,Y 



PQO 



— hi 



.1/4 



1 GeV^ 

(71) 

consistently with the fact that Ti r h > ^pqo and, thus, we 
obtain MD followed by RD during the era between the 
end of FHI and the onset of PQPT. On the other hand, 
Nhi* can be found from 



A ffl * 



A, 



HI 



A* 



(72) 



where Ahi is the total number of e-foldings generated 
during FHI and A* is the number of e-foldings elapsed 
from the onset of FHI until the scale fc* crosses outside 
the horizon of FHI and corresponds to the field value cr* . 
b. The power spectrum P-r. of the curvature per- 
turbation, which is generated during FHI and can be cal- 
culated at the pivot scale fc* as a function of cr* , is to be 
confronted with the WMAP7 data [l7J], i.e. 



pl/2 _ 



1 



K 



3/2 



HI 



2a/35 



|Vkl,<r| 



4.93 x 10 



-5 



(73) 



c. The (scalar) spectral index n s , which is given by 

n s = 1 - 6e(cr*) + 2*7(0-*), (74) 

is to be consisted with the fitting of the WMAP7 data 
by the standard power-low ACDM: 



n s = 0.963 ±0.028 



0.935 < ti s < 0.991 



(75) 



at 95% confidence level (c.L). Note, in passing, that the 
running of n s and the scalar-to-tcnsor ratio turn out to 
be vanishingly small in our model as in any model of FHI. 
d. In order for the PQPT to take place after a short 
temporary domination of Vpqo, the value of \P\ at the 
onset of PQPT, |-PpQi|, must be adequately larger than 
its critical value along the F-flat direction in Eq. (fTT|) . 
This entails 



|-PpQi| > /a/2 => s PQi > f a /Ap, 



(76) 



which implies mainly a lower bound on the parameter A 
through Eqs. (f6"2"|) and (|6"4")l . As we emphasize in Sec. Ill Al 
this requirement is directly related to the reduction of 
the G yield at the onset of nucleosynthesis, Y~, to an 
acceptable level. Indeed, Y~ in the case that Vpqo does 
not dominate can be estimated [27j as 



Y~ ~ 1.9 x 10~ 12 (Ti rh /10 10 GeV) 



(77) 



However, in the opposite case, if we take into account 
the entropy produced after PQPT during the subsequent 
reheating process - for computational details see Ap- 
pendix 0- we obtain 



K 



2G 



30 



Slrhi 



1/4 



T, 



2rh 



r T 



Y ~ 

1/4 1G 
PQO 



(78) 



which is suppressed relative to Y ~. In order to avoid 
spoiling the success of the SBB nucleosynthesis, an upper 
bound on Y~ is to be imposed depending on the G mass, 
to-, and the dominant G decay mode. For the conserva- 

G 

tive case that G decays with a tiny hadronic branching 
ratio, we have 



Y~< 

G ~ 



1Q -14 

10" 13 for 77i- 

io- 12 



0.69 TeV 

10.6 TeV (79) 

13.5 TeV 



respectively. The bound above can be somehow relaxed 
in the case of a stable G. However, it is achievable in our 
model, as we see below. 



VI. NUMERICAL RESULTS 

As can be easily seen from the relevant expressions 
above, our cosmological set-up depends on the following 
parameters: 

K, n a , A, M, f ai A M , n, a, b, c, d, e, /, and g. 

We fix f a = 10 12 GeV and A M = 0.01 so as to obtain 
/i ~ 1 TeV. We also set ?i = 5. These three parameters 
are involved in the determination of Ti r h and T2rh ~ via 
Eq. (f6"2"|) and (|70[) respectively - and play no crucial role 
in the inflationary dynamics. We also choose K a = 0.01 
and A = 0.1. Variation of n a leads to a variation of s 
according to Eq. ([54")) without creating drastic changes 
in the inflationary predictions for M, A^hi*, and ti s . On 
the other hand, A controls crucially s and Ti r h through 
Eqs. (15~4"1) and (f6~4")l and the upper bound on the condition 
of Eq. (110b ). Finally, we fix throughout our numerical 
computation c = d = e = f = g = 0.1 since these pa- 
rameters - contrary to a and b - have a minor impact on 



the calculation of 



and 



As we show below, the 



selected values above give us a wide and natural allowed 
region of the remaining fundamental inflationary param- 
eters (k,M, a, and b). Besides the parameters above, in 
our computation, we use as input parameters the quan- 
tities A* and fun — I with f = a, s, and q. We set 
/ ~ (1.5-4.5) xlO 17 GeV so as to obtain N m ~ 70-100. 
We then restrict M and A* so that Eqs. (fTTj) and ([75]) are 
fulfilled. We check also if the violation of Eqs. (fTUb ) and 
Eq. (JTUb) occurs according to the desired order so that 
our scenario is realized successfully. Using Eqs. (|64|) and 
(|74p , we can extract spqi and n s and compare them with 
the requirements of Eqs. ([76"]) and ([75]) . For the solutions 
presented below, Eq. ([76]) is safely fulfilled. 



10 



> 
a 



6 



JlljlljlllllinillllKHIliilm,, 



10 r 



1 r 



0.1 



0.01 r 



xxxxxxxxxxxxxxxxx 



xxxxxxxxxxxxx> 



20 40 60 80 100 120 140 160 180 

N 



30 








25 


^ 




I : 


> 20 
«15 


i 


' 


i : 


2 10 
" 5 



-5 

( 


1 








I 


) 5 10 15 20 


25 30 






a(10 16 GeV) 





FIG. 1: The evolution of a (black lines), s (gray lines), 
and q (light gray lines) as evaluated via the solution of 
Eq. ((48)) as functions of N for k — 0.002, K a = —b = 
0.01, a = -0.0125, \ = c = d = e = f = g = 0.1, 
and cthh = shh = <?hh = 3 x 10 17 GeV (solid lines) or 
o"hh = shh = Qhh = 2.5 x 10 17 GeV (dashed lines). Crosses 
are obtained by applying our analytical approach for aim = 
SHii = qm = 3 x 10 17 GeV. 



The selected / is to be large enough so that s reaches 
the attractor of Eq. (|M|) . Under this assumption, our 
results are independent of the precise /, as can be clearly 
deduced from Fig. [TJ where we plot a (black lines and 
crosses), s (gray lines and crosses), and q (light gray lines) 
as functions of N for k = 0.002, a = -0.0125, b = -0.01, 
and cthh = she = She = 2.5 x 10 17 GeV (dashed 
lines) or ernii = she = Qhu = 3 x 10 17 GeV (solid 
lines and crosses). The lines are drawn by solving nu- 
merically Eq. (|48p. whereas crosses are obtained by em- 
ploying Eqs. ([53)1 and (15Tj)) . For both choices of ernii = 
She = QBii'Sj we obtain m 2 _ = —0.0136, m 2 + = 1.82, 
N m * = 52.5, M = 2.86 x 10 16 GeV, n s = 0.963, 
ami = 4.08 x 10 16 GeV, and s H if = 1.3 x 10 14 GeV al- 
though in the hrst [second] case we obtain N{ — 87.1 
[JVf = 175.5], where N{ is the number of e-foldings 
elapsed from the commencement of FHI until Eq. (|10a ) 
is violated. It is impressive that M can take a value ex- 
actly equal to the SUSY GUT breaking scale, Mqut = 
2.86 x 10 16 GeV, contrary to all other realizations of the 
standard FHI - cf. Refs. @, OJ, HI HI- Despite the 
fact that T irh = 4.7 x 10 13 GeV, s P Qi/A P f a ~ 23 and, 
therefore, a second reheating process is possible which 
results to T2 r h = 3 x 10 4 GeV. It is worth emphasizing 
that K~~9x 10~ 9 , whereas Y~ ~ 1.9 x 10~ 14 which 
is consistent with the constraint of Eq. (|7^|) . If we set 
A M = 0.05, we take fi ~ 5 TeV, T 2rh = 1.5 x 10 5 GeV, 
and Y ~ ~ 9.5 x 10 -14 which again falls into the ranges of 
Eq. (ff 7 !)) . Finally, we observe that our analytical findings 
on a are very close to the numerical ones for all iV's. On 
the contrary, analytical and numerical results on s con- 
verge mainly for N ~ N{. Let us clarify that the results 
presented in the following are derived exclusively by our 



FIG. 2: The evolution of a and s as evaluated via the solution 
of Eq. igHJ in the a — s plane for k = 0.002, n a = —b = 
0.01, a = -0.0125, A = c = d = e = / = 5 = 0.1, and 
o"Hn = Shh = 3 x 10 17 GeV (solid line) or cthh = shh = 
2.5 x 10 17 GeV (dashed line). The gray [light gray] region is 
excluded by Eq. flOk) [Eq. (fTUb)]. 



numerical program. 

As observed in Fig. [TJ immediately after the onset of 
FHI, q decreases sharply and so it does not influence the 
inflationary dynamics. The evolution of the other two 
fields (a and s) can be represented in the a — s plane 
as in Fig. [2) where we show the evolution of a and s for 
the parameters used in Fig. [1] and taking <7he = she = 
3 x 10 17 GeV (solid line) or cr H ii = Shu = 2.5 x 10 17 GeV 
(dashed line). There, we draw also the gray [light gray] 
region which is excluded by Eq. (|10a ) [Eq. (|10b )]. We 
remark that Eq. (fl"0a) is violated earlier and so the pre- 
ferred hierarchy in the domination of Uhio or Vpqo in our 
cosmological proposal is valid - see Sec. Ill Al 

The importance of the coefficient a in reducing n s can 
be easily concluded from Fig. Efa), where we depict n s 
versus a for b = —0.01 and various k's indicated in the 
graph. Increasing the absolute value of a, \m 2 _\ increases 
too and so n s decreases and becomes consistent with the 
observationally favored range of Eq. (|75l) . On the other 
hand, the same variation of a leads to a reduction of M 
which lies around its SUSY GUT value, Mqut, as can 
be concluded from Fig. [3Jb). 

Confronting FHI with all the constraints of Sec. |V] we 
can delineate the allowed (lightly gray shaded) region in 
the k — a [k — M] plane as in Fig. 0|a) [Fig. HJb)] , where 
we take b = —0.01 and show the adopted conventions 
for the various lines. In particular, the gray dashed [dot- 
ted] lines correspond to n s = 0.991 [n s = 0.935], whereas 
the gray solid lines are obtained by fixing n s — 0.963 
- see Eq. ((75)) . The black solid lines correspond to 
M = Mqut, whereas beyond the black dashed lines 
the preferred hierarchy in the violation of Eqs. ((TOr ) and 
((TUb ) fails. We observe that the latter requirement holds 
only for n < A. Along the gray solid line (n s = 0.963), 
we obtain 2.7 x 10 13 < Ti rh /GeV < 1.5 x 10 14 and 
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and several «'s indicated in the graphs. The region of Eq. Q75p is also limited by thin lines in (a). 
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1.9 x 1(T 14 < Y~ < 1CT 13 for 2 x 10" 4 < k < 0.113. 
Note that T2 r h = 3 x 10 4 GeV remains fixed since A M and 
f a are also fixed throughout our computation. 

One of the outstanding features of our proposal is that 
the reduction of n s can be attained without disturbing 
the monotonicity of the potential, contrary to other sim- 
ilar suggestions - see Refs. [20, [23, |24 |. This fact can 
be highlighted in Fig. [SJ where we present the variation 
of the inflationary potential Vhi as a function of the in- 
flaton field, a, for the values of the parameters corre- 
sponding to the three intersections of the solid black line 
with the gray lines in Fig.[D[a) and (b). Namely, we take 
M = Mgut, b = -0.01, and k = 0.0036, a = -0.0062, 
and CTHii = s H ii = gmi = 3 x 10 17 GeV (n s = 0.991, 
dashed line) or k = 0.002, a = —0.0125, and cthh = 
she = SHii = 2.5 x 10 17 GeV (n s = 0.963, solid line) or 
k = 0.00111, a = —0.0189, and ohe = she = <7Hii = 



2 x 10 17 GeV (n s = 0.935, dotted line). The values cor- 
responding to o* and Ui are also designed. We observe 
that for large values of a, Vhi develops an oscillatory 
behavior due to the initial oscillations of s and q - see 
Fig. [1] However, Vm, for lower ct's, remains monotonic 
and, therefore, no complications arise in the realization 
of the inflationary dynamics. 

Letting b vary for a number of fixed values of a and for 
the most exciting case with M = Mgut, we can depict 
the values allowed by all the constraints of Sec. |V]in the 
k — b plane - see Fig. [5] The various lines terminate 
at low [high] n's due to the saturation of Eq. ([75"|) from 
below [above]. The central n s is obtained at n ~ 0.002. 
We readily conclude that the allowed a's and b's for M — 
Mqut and fixed n s are almost K-independent. This is 
because, for fixed n s , m_ is fixed too. In particular, for 



0.963, we have m_ 



-0.0136, whereas, for n s 
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FIG. 5: The variation of the inflationary potential Vhi as a 
function of a for M = Mgut, K a — — b — 0.01, A = c = d — 
e = / = 9 = 0.1, and re = 0.0036, a = -0.0062, and ct H h = 
SHli = Qhh = 3 x 10 17 GeV (n s = 0.991, dashed line) or re = 
0.002, a = -0.0125, and cthh = Shh = «hh = 2.5 x 10 17 GeV 
(n s = 0.963, solid line) or k = 0.00111, a = -0.0189, and 
cthii = Shu = qmi = 2 x 10 17 GeV (n s = 0.935, dotted line). 
The values corresponding to <r, and erf are also depicted. 



0.935 [n s = 0.991], we have ml ~ -0.02 [ml ~ -0.008] 
and « ~ 0.0036 [k ~ 0.0011]. In all cases, m\ ~ 1.82. 
For this reason, we can present in Fig. [7] the allowed 
values by Eqs. ([7T]>. l(75|) . and ((75J in the a - 6 plane 
for M ~ Mgut and k ~ 0.002 (n s = 0.963, solid lines), 
K ~ 0.0036 (ns = 0.991, dashed lines) or k ~ 0.0011 
(n s = 0.935, dotted lines). We observe that our scenario 
can be realized for both signs of a and 6, contrary to 
the cases studied in Refs. [20, [23, [24[ where negative 6's 
are necessitated. Moreover, compared to the latter cases, 
larger |6|'s (of the order of 0.1) are here permitted. 
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FIG. 6: Allowed values by Eqs. |[7I)|. (75) , and (J75J) in the 
k — b plane for K a = 0.01, X = c = d — e — f — g — 0.1, M = 
A^gut, and various a's indicated on the curves. 



FIG. 7: Allowed values by Eqs. (JTTJl , ([73), and J75) in the 
a — & plane for K a = 0.01, X — c — d — e — f — g — 0.1, M ~ 
A/gut, and n — 0.002 (solid lines), re ~ 0.0036 (dashed lines) 
or re ~ 0.0011 (dotted lines). 



VII. CONCLUSIONS 

We investigated a cosmological scenario according to 
which a SUSY GUT scale FHI is followed by a PQPT 
which resolves the strong CP and the \i problems of 
MSSM. The PQPT is tied to renormalizable superpo- 
tential terms and the possible catastrophic production of 
domain walls can be eluded by the introduction of extra 
matter superfields which can be chosen so that the MSSM 
gauge coupling constant unification is not disturbed. The 
inflaton-like field, associated with PQPT, plays a crucial 
role in the construction of the Kahler potential which 
is expanded up to fourth order in powers of the various 
fields. The FHI reproduces the current data on P-r and 
n s within the power-law ACDM cosmological model and 
generates the number of e-folds required from the reso- 
lution of the horizon and flatness problems of the SBB. 
The dynamics of FHI is investigated both numerically 
and analytically and the results are compared with each 
other. Fixing n s to its central value and M to the SUSY 
GUT scale, we concluded that k — 0.002 with the re- 
maining parameters taking more or less natural values, 
±(0.01 — 0.1). It is gratifying that our model supports a 
second stage of reheating after PQPT, which dilutes suffi- 
ciently the G abundance so as to become observationally 
safe for G masses even lower than 10 TeV. 
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Appendix A: Reheating Processes and Gravitino 
Constraint 

In this Appendix, we present a numerical description 
of the post-inflationary evolution of the various energy 
densities in our set-up paying special attention to the 
dilution of the G yield, Y~. 

The energy density, p\ [p 2 ], of the oscillatory system 
which reheats the universe at the temperature Ti r h p2 r h] , 
the energy density of produced radiation, pp,, and the 
number density of G, np^, satisfy the following Boltzmann 
equations - cf. Refs. [23ll27|: 



Pi- 
f>2- 
Pr 

it- 



3H Pl +T lPl = Q, 
3Hp 2 + T2P2 = 0, 
- 4Fp R - Tip! - r 2 p 2 = 0, 



3Hn- 



C^in^f 



0. 



(Al) 
(A2) 
(A3) 
(A4) 



Here n cq = C(3)T 3 /-7r 2 is the equilibrium number den- 
sity of the bosonic relativistic species, C~ is a collision 

term for G production whi ch, in the limit of the massless 
gauginos, turns out to be [23, [5l| 



C 



3tt 



G 



16C(3)m 2 ^ 



3 



c l g l In 



(A5) 



where gi (with i = 1,2,3) are the gauge coupling con- 
stants of MSSM, (ki) = (1.634, 1.312, 1.271), and (a) = 
(33/5, 27, 72). Also the Hubble expansion parameter, H, 
during this period is given by 



H 



1 



n/3: 



mp 



m~n~ 



Pl+ P2+ PR 



1/2 



(A6) 



Clearly, in the limit of massless MSSM gauginos, the n~ 
computation is m —-independent. The temperature, T, 
and the entropy density, s (not to be confused with the 
field s), can be found using the relations 



PR 



:-!() 



<7*T and s 



45 



■g*T\ 



(A7) 



where g*(T) = gi r h* [g*(T) = #2rii*] for T > T PQ 
[T < Tpq] with Tpq being defined as the solution of the 
equation pr (Tpq) — Vpqo an d can be found numerically. 
The numerical integration of Eqs. (|A1I) — (|A4|) is facili- 
tated by absorbing the dilution terms. To this end, we 
find it convenient to define [47j the following dimension- 
less variables 



h = piR\ h = P2R 3 , /r = prR 4 



and f~ = n~R 
■' G G 



Converting the time derivatives to derivatives w.r.t N = 
In (R/Rnif ) with i?nif being the value of the scale factor 
at the end of FHI (the value of i?mf turns out to be 
numerically irrelevant), Eqs. (|Alj) - (|A4j) become 



Hf[ 



-ri/i 



(A9) 
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line), log Vpqo (black dashed line), and logY— (black solid 
line) as functions of logT for k = 0.002, n a — —b — 0.01, a = 
-0.0125, A = c = d = e = / = 3 = 0.1, M = M GU t, fa = 
10 12 GeV, and A M = 0.01. 



Hf 2 -- 


= -r 2 / 2l 


(A10) 


Hfk - 


= TJ.R + TzhR, 


(All) 


Hf~ - 


= C~(n^fR\ 


(A12) 



where the differentiation w.r.t. N and N - see Sec. IIIIBI 
- coincides. Also H and T can be expressed in terms of 
the variables in Eq. (|A8|) as 



H 



i 5 f S + h+f2 + fn/R 



andT 



30 / R 



ir 2 g*R 4 
(A13) 
The system of Eqs. (|A9|) — (IA12|) can be solved, imposing 
the following initial conditions (the quantities below are 
considered functions of the independent variable N): 

Pi(0) = Vbo, p R (0) - n 5 (0) = 0, and p 2 (A PQ ) = V PQ0 , 

where A^pq is the value of N corresponding to the tem- 
perature Tpq. Needless to say that we set P2(N) = for 

^V < TVpQ. 

In Fig. [5J we illustrate the cosmological evolution of 
the quantities log p.; with i = 1 (dotted gray line), i = 2 
(dashed gray line), and i — R (gray line), log Vpqo 
(black dashed line), and logY- (black solid line) as func- 
tions of log T for the values of the parameters adopted 
in Fig. [T] In particular, the parameters which deter- 
mine the evolution of the various quantities during the 
(A8) post-inflationary period are K,X,M,K a ,f a , and A M . We 



take k = 0.002, A = 0.1, M = M GVTl n a = 0.01, f a = 
10 12 GeV, and A M = 0.01. From Fig. H we observe that 
FHI is followed successively by the following four epochs: 
(i) a MD era, due to the oscillating and decaying inflaton 
system, which lasts until T ~ Ti r h given by Eq. (l6"2")) , (ii) 
a RD epoch, which terminates at Tpq, (iii) a MD era ere- 
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ated by the oscillations of the PQ-system, which is com- 
pleted at T ~ T 2r h given by Eq. ([70]) , and (iv) a RD epoch 
after which the universe enters the SBB phase. The com- 
pletion of the two reheating processes corresponds to the 
two intersections of p\ and p 2 with pp i n Fig. [8j Here 
we omit for simplicity any possible instantaneous domi- 
nation of Vpqo over the various energy densities since, as 
we mention in Sec. IIV1 Vpq does not support any signif- 
icant period of inflation. In Fig. [8j we also see that the 
G yield, Y~ = n~/s, takes the value Y ± ~ ~ 1.1 x 1CT 8 

- see Eq. (|77j) - for T ~ Ti r h- However, due to the en- 
tropy released during the out-of-equilibrium decay of the 
PQ system, Y~ decreases sharply to Y ~ ~ 3.3 x 10~ 14 

- see Eq. (j78|) . This result can be understood from the 



following relation 



y _y 5(Tp Q ) /i?p Q \ : 

-'■ 1G s(T 2rh ) \R 2lh J 



2G 



(A15) 



where Rpq [i?2rh] is the value of the scale factor cor- 
responding to TpQ [T 2r h]. Taking into account Eq. (|A7[) 

— 1/3 

and the fact that, during this MD era, R oc p osc - recall 
that /o osc is the energy density of the oscillating system 
- we arrive at Eq. (1781) . We verify that the analytical 
results agree remarkably with the numerical ones. 

Therefore, we can easily appreciate the importance of 
PQPT, as realized in our scenario, in lowering Y~ to a 

1 Ct 



value compatible with the observational data 27 1 
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